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Abstract 

Background: Although 78X7 mutations have been identified in patients with 22q11.2 deletion syndrome (22q1 1.2DS)-like 
phenotypes including characteristic craniofacial features, cardiovascular anomalies, hypoparathyroidism, and thymic 
hypoplasia, the frequency of 78X7 mutations remains rare in deletion-negative patients. Thus, it would be reasonable to 
perform a comprehensive genetic analysis in deletion-negative patients with 22q1 1.2DS-like phenotypes. 

Methodology/Principal Findings: We studied three subjects with craniofacial features and hypocalcemia (group 1), two 
subjects with craniofacial features alone (group 2), and three subjects with normal phenotype within a single Japanese 
family. Fluorescence in situ hybridization analysis excluded chromosome 22q11.2 deletion, and genomewide array 
comparative genomic hybridization analysis revealed no copy number change specific to group 1 or groups 1+2. However, 
exome sequencing identified a heterozygous 78X7 frameshift mutation (c.1253delA, p.Y418fsX459) specific to groups 1+2, 
as well as six missense variants and two in-frame microdeletions specific to groups 1+2 and two missense variants specific to 
group 1. The 78X7 mutation resided at exon 9C and was predicted to produce a non-functional truncated protein missing 
the nuclear localization signal and most of the transactivation domain. 

Conclusions/Significance: Clinical features in groups 1+2 are well explained by the 76X7 mutation, while the clinical effects 
of the remaining variants are largely unknown. Thus, the results exemplify the usefulness of exome sequencing in the 
identification of disease-causing mutations in familial disorders. Furthermore, the results, in conjunction with the previous 
data, imply that 78X7 isoform C is the biologically essential variant and that 78X7 mutations are associated with a wide 
phenotypic spectrum, including most of 22q11.2DS phenotypes. 
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Introduction 

Chromosome 22qll.2 deletion syndrome (22qll.2DS) is a 
developmental disorder associated with characteristic craniofacial 
features with velopharyngeal incompetence, cardiovascular anom- 
alies primarily affecting the outflow tracts, hypoparathyroidism 
and resultant hypocalcemia, and thymic hypoplasia leading to 
susceptibility to infection [1]. This condition is also frequently 
accompanied by non-specific clinical features such as develop- 
mental retardation [1]. Expressivity and penetrance of these 
features are highly variable and, consistent with this, chromosome 
22qll.2 deletions have been identified in DiGeorge syndrome 



(DGS) and velocardiofacial syndrome (VCFS) with overlapping 
but different patterns of clinical features [1] . 

While multiple genes are involved in chromosome 22qll.2 
deletions [2], TBX1 (T-box 1) has been regarded as the major gene 
relevant to the development of clinical features in 22ql 1.2DS [3]. 
Indeed, heterozygous TBX1 mutations have been identified in 
several deletion-negative patients with 22qll.2DS phenotype [2- 
8], and mouse studies argue for the critical role of Tbxl in the 
development of 22qll.2DS phenotypes [3]. However, the 
frequency of TBX1 mutations remains rare in deletion-negative 
patients: Gong et al. identified only a few probable TBX1 
mutations after studying 40 patients with DGS/VCFS phenotypes 
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[4], and Zweier et al. found a single TBX1 mutation after 
examining 10 patients with 22qll.2DS phenotype [8]. This 
indicates the presence of genetic heterogeneity in the development 
of 22qll.2DS phenotype in deletion-negative patients. Consistent 
with this, another DGS/VCFS locus has been assigned to 
chromosome 1 Op 13- 14 region [9]. Thus, it would be reasonable 
to perform a comprehensive genetic analysis in deletion-negative 
patients with 22qll.2DS phenotype. 

In this regard, recent advance in molecular technologies has 
enabled to perform comprehensive genetic analyses, thereby 
contributing to the identification of underlying factors in genetic 
disorders. Indeed, genomewide array comparative genomic 
hybridization (CGH) has identified multiple disease-associated 
copy- number changes [10], and exome sequencing has discovered 
multiple disease-causing gene mutations [11]. In particular, these 
technologies can be powerful methods for familial disorders, 
because it is predicted that a single copy-number change or 
mutation is shared in common by affected subjects and is absent 
from non-affected subjects within a family. 

Here, we performed array CGH analysis and exome sequencing 
in a family with 22ql 1.2DS-like clinical features. Although this 
study did not discover a novel disease gene, a TBX1 mutation was 
successfully identified. 

Materials and Methods 

Ethics statement 

The Institutional Review Board Committees of Hamamatsu 
University School of Medicine, Tohoku University School of 
Medicine, Kurashiki Central Hospital, and National Research 
Institute for Child Health and Development considered and 
approved the study, consent/assent procedures, and the publica- 
tion of images and case details associated with this work. The 
individuals in this manuscript have given written informed consent 
(as outlined in PLOS consent form) to publish these case details. 
Actually, this study was performed after obtaining written 
informed consent from the parents of the child subjects and from 
the adult subjects. Furthermore, the mother and the elder brother 
aged 19 years old have given written informed consent to 
publication of the facial photographs of the two brothers; in 
addition, the younger brother aged 10 years has given informed 
assent. 

Clinical Report 

The pedigree of this Japanese family is shown in Fig. 1 , and 
clinical findings of the family members are summarized in Table 1 . 
The proband (subject III-5) was found to have hypocalcemia and 
hyperphosphatemia in a pre-operation laboratory test for repeated 
otitis media at 8 years of age, and was referred to Department of 
Pediatrics at Kurashiki Central Hospital. Subsequent examination 
revealed borderline low serum intact PTH value. Thus, he was 
diagnosed as having hypoparathyroidism, and received vitamin D 
therapy. Furthermore, physical examination showed characteristic 
craniofacial features with velopharyngeal incompetence suggestive 
of22qll.2DS. 

Similar craniofacial features were also exhibited by subjects II-2, 
III- 1 , III-6, and III-7, and hypocalcemia was also identified in 
subjects II-2 and III-7. Actually, subject II- 2 was taking vitamin D, 
and subject III-7 was noticed to have hypocalcemia at birth 
because of the history of subject III-5, and was treated with 
vitamin D. The five subjects with 22ql 1.2DS-like craniofacial 
features lacked cardiovascular anomalies; while they also lacked 
susceptibility to infection, except for repeated otitis media in 
subject III-5, thymic hypoplasia was not evaluated in four of the 



five subjects. By contrast, the five subjects exhibited borderline to 
mild developmental delay. Indeed, adult subjects II- 2 and III- 1 
had some difficulty in verbal communications, although they were 
able to get on their daily life, and subject II-2 was able to take care 
of family members. Similarly, child subjects III-5, III-6, and III-7 
also showed speech delay, and subjects III-5 and III- 7 received 
speech therapy. Furthermore, subject III-7 was diagnosed as 
having pervasive developmental disorder, and his verbal, perfor- 
mance, and full scale intelligence quotients were assessed as 63, 64, 
and 60, respectively, by the WISC-III method at 10 years of age. 
In addition, subject II-2 had sensorineural deafness, and subject 
III- 5 had Graves' disease. 

Molecularly Studied Subjects 

Molecular studies were performed for eight subjects in this 
family, using peripheral blood samples. They were divided into 
three groups in terms of clinical findings: group 1, subjects II-2, 
III-5, and III-7 with craniofacial features and hypocalcemia; group 
2, subjects III- 1 and III-6 with craniofacial features alone; and 
group 3, subjects II- 1, III- 3, and IV- 1 with apparendy normal 
phenotype (Fig. 1). 

FISH and Array CGH Analyses 

Fluorescence in situ hybridization (FISH) analysis was performed 
with a probe for HIRA on the commonly deleted chromosome 
22ql 1.2 region and that for ARSA at chromosome 22ql3 utilized 
as an internal control (Abott). Array CGH was carried out using a 
genomewide 2x400K Agilent platform catalog array, according to 
the manufacturer's instructions (Agilent Technologies), and copy 
number variants/polymorphisms were screened with Agilent 
Genomic Workbench software using the Database of Genomic 
Variants (http:/ /dgv.tcag.ca/ dgv/ app/home). 

Exome and Sanger Sequencings 

Exon capture was performed with the SureSelect Human All 
Exon kit v4 (Agilent Technologies). Exon libraries were sequenced 
with the Illumina Hiseq 2000 platform according to the 
manufacturer's instructions (Illumina), providing 108-122 average 
depth for each sample. Paired 101 -base pair reads were aligned to 
the reference human genome (UCSChgl9) using the Burrows- 
Wheeler Alignment tool [12]. Likely PCR duplicates were 
removed with the Picard program (http://picard.sourceforge. 
net/). Single-nucleotide variants and indels were identified using 
the Genome Analysis Tool Kit (GATK) vl.6 software [13]. SNVs 
and indels were annotated against the RefSeq database, 1000 
Genomes Project variant data, and dbSNP135 with the ANNO- 
VAR program [14]. 

To confirm mutations indicated by exome sequencing, Sanger 
sequencing was performed for PCR products obtained with 
primers flanking the detected mutations, using a 3500xL genetic 
analyzer (Life Technologies). Furthermore, the PCR products 
were subcloned with TOPO TA Cloning Kit (Life Technologies), 
and normal and mutant alleles were sequenced separately. 

In silico protein functional analysis 

Function of proteins with missense variations was assessed by 
Polymorphism Phenotyping-2 (PolyPhen-2, http://genetics.bwh. 
harvard.edu/pph2/) and Sorting Intolerant From Tolerant (SIFT, 
http://sift.jcvi.org/), and that of proteins with in-frame amino 
acid deletions was evaluated by PROVEAN predictions (http:// 
provean.jcvi.org/index.php). 
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Figure 1. The pedigree of this family. Facial features of subjects 111-5 and 
doi:10.1371/journal.pone.0091598.g001 



111-7 

(3.7 years) 

1-7 are shown. 



Results 

FISH and Array CGH Analyses 

FISH analysis delineated two signals for HIRA (Fig. 2A). Array 
CGH analysis revealed no copy number change specific to group 1 
or groups 1+2, in the entire genome including chromosome 
1 Op 13-14 and chromosome 22ql 1.2 regions (Fig. 2B). 

Exome and Sanger Sequencings 

Exome sequencing identified nine heterozygous non-synony- 
mous variants (six missense variants, two in-frame microdeletions, 
and one frameshift variant) that were specific to groups 1+2 
(namely, they were present in groups 1+2 and absent from group 3 
as well as from 1000 Genomes, dbSNP135, and our in-house 
exome data from 70 individuals) (Table SI). Notably, the 
frameshift variant (c.l253delA, p.Y418fsX459) was found at exon 
9C of TBX1 for DGS/VCFS (Fig. 3). Of the remaining eight 
variants, two variants were also detected in disease-related genes: 
p.G204R in HDAC4 for brachydactyly-mental retardation syn- 
drome [15], and p.276del in CCND1 constituting a susceptibility 
factor for colorectal cancer and a modifier for von Hippel-Lindau 
disease [16,17]. Exome sequencing also detected two heterozygous 
missense variants that were specific to group 1 (Table SI). 

When all variants were included, exome sequencing revealed: 
(1) 83 non-synonymous and 86 synonymous variants that were 
present in groups 1+2 and absent from group 3 (Table S2); (2) 54 
non-synonymous and 48 synonymous variants that were present in 
group 1 and absent from groups 2+3 (Table S3); (3) 6,033 non- 
synonymous and 6,667 synonymous variants that were present in 
groups 1+2, but not specific to groups 1+2 (thus, they may be 
present in group 3 or absent from group 3); and (4) 7,073 non- 
synonymous and 7,861 synonymous variants that were present in 
group 1 , but not specific to group 1 . Furthermore, comparison of 
the exome sequencing data between group 1 with hypocalcemia 
and group 2 without hypocalcemia revealed 23 1 non-synonymous 
and 254 synonymous variants that were present in group 1 and 
absent from group 2 (Table S4), and 246 non-synonymous and 
242 synonymous variants that were present in group 2 and absent 
from group 1 (Table S5). (The variant data other than those 
described in Supplemental Tables may be available on request 



after discussion with the family members and approval by our 
IRBs, because they contain a huge amount of individual genetic 
information.) 

In silico protein functional analysis 

The results are summarized in Table SI. The p.G204R in 
HDAC4 and the p.E276del in CCMD1 were assessed as non- 
pathologic, while some variants were evaluated as potentially 
pathologic. 

Discussion 

Exome sequencing successfully identified a heterozygous 
frameshift variant on exon 9C of TBX1. The c.!253delA 
(p.Y418fsX459) appears to be a disease-causing mutation, because 
it is predicted that this variant escapes nonsense-mediated mRN A 
decay due to its position on the final exon [18] and produces a 
truncated protein lacking the nuclear localization signal (NLS) and 
most of the transactivation domain (TAD) on exon 9C (Fig. 3) 
[19]. In support of this, functional studies for a similar c.l223delC 
(p.S408fsX459) mutation on exon 9C have shown that the 
truncated p.S408fsX459 protein was incapable of localizing to 
nucleus and lost transactivation function [2,5,19]. One may argue 
that this c.l253delA mutation affects TBX1 isoform C (TBX1C, 
TBX 1-003) alone, while TBX1 produces three transcript variants 
containing T-box [2,4] (Fig. 3). However, TBX1C is the major 
transcript with the NLS and the TAD in human and is highly 
homologous to mouse Tbxl [4] (Fig. SI). 

Craniofacial features in groups 1+2 and hypocalcemia in group 
1 are well explained by the TBX1 mutation [3] . This argues for a 
critical role of this mutation in the phenotypic development in 
groups 1+2, while the clinical effects of the remaining variants 
identified by exome sequencing are largely unknown. In this 
regard, comparison between group 1 with hypocalcemia and 
group 2 without hypocalcemia revealed a large number of non- 
synonymous and synonymous variants that exclusively belonged to 
either group 1 (Table S4) or group 2 (Table S5), although the lists 
did not contain a c.2968A>G (p.R990G) SNP in CASR (calcium 
sensing receptor) that has a gain-of-function effect and appears to 
raise the susceptibility to hypocalcemia (Fig. S2) [20]. Thus, it is 
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Table 1. Clinical findings of the family members. 







TBX1 r 


nutation (+) 








TBX1 1 


nutation (— ) 




TBX1 r 


nutation (N.E. 


) 


Individual 


11-2 


111-1 


■11-5 


II 1-6 


111-7 


11-1 


■ 11-3 


IV-1 


11-3 


11-4 


111-2 


Present age (year) 


51 


26 


19 


13 


10 


59 


22 


5 


50 


49 


25 


Sex 


F 


F 


M 


F 


M 


M 


F 


M 


M 


M 


M 


Craniofacial features 


+ 


+ 


+ 


+ 


+ 














Hypertelorism 


+ 


+ 


+ 


+ 


+ 














Blepharophimosis 


+ 


+ 


+ 


+ 


+ 














Low set ears 


+ 


+ 


+ 


+ 


+ 














Auricular anomalies 


+ 






















Narrow nose 


+ 


+ 


+ 


+ 


+ 














Cleft palate 
























Micrognathia 


+ 


+ 


+ 


+ 


+ 














Velopharyngeal incompetence 




+ 


+ 


+ 


+ 














Hypoparathyroidism 


+ 




+ 




+ 














Age at examination (year) 


44 


17 


8 


4 


0 (1 day) 


N.E. 


15 


0 (6 days) 


N.E. 


N.E. 


18 


Serum calcium (mg/dl_) a 


7.6 e 


9.0 


6.0 


9.1 


5.9 




9.0 


9.8 






9.6 


Serum i-phosphate (mg/dl_) a 


3.9" 


4.9 


9.1 


5.0 


N.E. 




4.8 


6.3 






4.6 


Serum intact PTH (pg/dl_) a 


31 e 


N.E. 


15 


N.E. 


19 




N.E. 


34 






N.E. 


Cardiovascular anomalies' 5 
























Hypoplastic thymus c 




N.E. 


N.E. 


N.E. 


N.E. 


N.E. 


N.E. 


N.E. 


N.E. 


N.E. 


N.E. 


Susceptible to infection 






_f 


















Other features 
























Developmental retardation 


+ 


+ 


+ 9 


+ 


+ 9 














Sensorineural deafness 


+ h 






















Graves' disease 






+' 



















Individuals correspond to those shown in Fig. 1. 

i-phosphate: inorganic phosphate; SD: standard deviation; F: female; M; male; and N.E.: not examined. 

a Reference values: calcium, 9.0-11.0 mg/dL in infants and 8.8-10.2 mg/dL in adults; inorganic phosphate, 4.8-7.5 mg/dL in infants and 2.5-4.5 mg/dL in adults, and 
intact PTH, 10-65 pg/dL in infants and 14-55 pg/dL in adults. 

Conversion factor to the SI unit: 0.25 for calcium (mmol/L), 0.32 for inorganic phosphate (mmol/L), and 0.106 for intact PTH (pmol/L). 

b Examined by echocardiography, chest roentgenography, and/or electrocardiography. 

c Examined by computed tomography. 

d Received velopharyngeal closure. 

e On treatment with vitamin D. 

'Repeated otitis media only. 

9 Received speech therapy. 

h Required hearing aids. 

'At the time of diagnosis (1 1 years of age), serum TSH was <0.01 mlU/L, free T 3 33.1 pg/mL [51.0 pmol/L], free T 4 5.1 1 ng/dL [65.8 nmol/L], and TSH receptor antibody 

1284% [normal range <1.9%]. 

doi:1 0.1 371 /journal.pone.0091 598.t001 



likely that, together with environmental factors, the combination 
of hitherto unknown calcium metabolism-related functional 
variants would underlie different serum calcium values between 
groups 1 and 2. 

In addition to craniofacial features with and without hypocal- 
cemia, TBX1 mutation positive subject II-2 had sensorineural 
deafness, and III-5 had Graves' disease. Since such features are 
occasionally manifested by patients with 22qll.2DS [21,22], the 
results may suggest the relevance of TBX1 to such rather 
infrequent features in 22qll.2DS. 

The five TBX1 mutation positive subjects in groups 1+2 lacked 
cardiovascular lesion and manifested borderline to mild develop- 
mental retardation (while they had no susceptibility to infection, 
assessment of thymic hypoplasia remained fragmentary). By 
contrast, cardiovascular lesion is frequently observed and devel- 
opmental retardation is rare in previously reported patients with 



TBX1 mutations, although clinical features are fairly variable 
among mutation positive patients (Table 2). Such difference would 
more or less be ascribed to an examination bias that TBX1 has 
been analyzed in patients with isolated cardiovascular lesion in 
several studies [4,6,7] or to the functional difference of the mutant 
proteins [2,5-8]. However, in seven patients who have been 
examined for DGS/VCFS-like clinical features and found to have 
frameshift mutations on exon 9C (p.S408fsX459, p.H425fsX613, 
and p.S431fsX608) affecting the NLS and the TAD, cardiovas- 
cular lesion was present in four patients and developmental delay 
was absent or not described, despite apparent similarity in the 
ascertainment of patients and the function of mutant proteins 
between the seven patients and the five affected subjects in this 
family (Table 2) [2-5]. 

Thus, there may be protective factor(s) for cardiovascular lesion 
and susceptibility factor(s) for developmental delay in groups 1+2. 
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-2 
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Figure 2. FISH and array CGH analyses in the proband ( 1 11-5). A. 

FISH analysis. Two signals are shown for both HIRA at 22q11.2 (red 
signals indicated by arrows) and ARSA at 22q13 (green signals indicated 
by arrowheads). B. Array CGH analysis. No copy number change is 
found for chromosome 10 carrying the second DiGeorge region and 
chromosome 22 harboring the DGS/VCFS critical region, as well as other 
chromosomes (not shown). Black, red, and green dots denote signals 
indicative of the normal, the increased (>+0.5), and the decreased 
(<— 0.8) copy numbers, respectively. Although several red and green 
signals are seen, there is no portion associated with >3 consecutive red 
or green signals. 

doi:1 0.1 371 /journal.pone.0091 598.g002 

In this regard, a simple explanation would be that protective 
factor(s) for cardiovascular lesion are present in groups 1+2 and 
may be present in group 3 or absent from group 3, whereas 
susceptibility factor(s) for developmental delay is present in groups 
1+2 and absent from group 3. Since 6,033 non-synonymous and 
6,667 synonymous variants were found to be present in groups 



1+2 but not specific to groups 1+2, and 83 non-synonymous and 
86 synonymous variants were revealed to be present in groups 1+2 
and absent from group 3, a certain fraction of functional variants 
may constitute protective factor(s) for cardiovascular lesion and 
susceptibility factor(s) for developmental delay. In addition, while 
p.G204R on HDAC4 for brachydactyly-mental retardation syn- 
drome was assessed as non-pathologic by in silico analysis, it may 
have played a certain role in the occurrence of developmental 
delay in groups 1+2. Actually, such protective and susceptibility 
factor(s) would be more complex, with the effects of functional 
variants unique to each patient as well as the influences of 
environmental factors. Furthermore, it remains possible that the 
c.l253delA (p.Y418fsX459) mutation found in this study may be 
related to a specific phenotype characterized by the presence of 
craniofacial features and developmental delay and by the absence 
of cardiovascular lesion, because of a hitherto unrevealed 
mechanism(s). This matter awaits further studies. 

Besides the clinical findings, several matters are also notable in 
the nine apparently pathologic TBX1 mutations identified to date 
(Table 2). First, the mutations reside on exons 3-8 common to 
isoforms A-C or on exon 9C specific to isoform C, with no 
mutation on exons 9A and 9B specific to isoforms A and B. This 
would be consistent with TBX1C having the primary biological 
function. Second, while most mutations have loss-of-function 
effects, gain-of-function effects have been suggested for p.F148Y, 
P.H194Q, and p.310S by in vitro studies [8]. Thus, TBX1 loss-of- 
function mutations and gain-of-function mutations may result in 
overlapping clinical features. Lastly, the c.l274_1281delAC- 
TATCTC (p.H425fsX613) missing the NLS on exon 9C was 
shared by a patient with DGS-like phenotype and the apparently 
normal mother, and the c.l29_185del57 (p.43-61dell9) with 
reduced transcriptional activity was common to a patient with 
non-syndromic tetralogy of Fallot and the apparently normal 
mother. This would imply the reduced penetrance of phenotypes 
caused by these mutations. 

In summary, we identified a TBX1 mutation by exome 
sequencing in a family with chromosome 22qll.2 deletion-like 
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c. 1 253delA (p. Y41 8fsX459) 



Figure 3. TBX1 mutation identified in this family. A. Genomic structure of TBX1 and the position of the mutation. The color and the white boxes 
represent the coding regions and the untranslated regions on exons 1-10 (El —El 0), respectively; the red, the purple, and the orange segments 
indicate the coding regions on the final exons 9C, 9A, and 9B (splice variants), respectively. The T4dox is indicated by yellow boxes, the nuclear 
localization signal (NLA) by a blue segment, and the transactivation domain (TAD) by a green arrow. The c.1253delA (p.Y41 8fsX459) identified in this 
family resides on exon 9C. B. Transcripts of TBX1. Three variants are formed by alternative splicing of the final exons 9C, 9A, and 9B. The c.1 253delA 
(p.Y418fsX459) mutation is predicted to yield a truncated TBX1C protein missing the NLS and most of the TAD. The stippled box of p.Y418fsX459 
denotes aberrant amino acid sequence produced by the frameshift mutation. C. Electrochromatograms showing the frameshift mutation by Sanger 
sequencing. The primer sequences used are: 5'-GCGGCCAAGAGCCTTCTCT-3' and 5'-GGGTGGTAGCCGTGGCCA-3'. 
doi:1 0.1 371 /journal.pone.0091 598.g003 
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Table 2. Summary of patients with TBX1 mutations. 
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Exon 3 
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c.1 223 


c.1253 


C.1274J281 


c.1293_1315 


c.1399_1428 


C.129J85 


c.443T>A 


C.5820G 


c.928G>A Deletion 




deIC 


delA 


del8 


del23 h 


dup30 k 


del57 k 








Amino acid 


p.S408 


p.Y418 


p.H425 


p.S431 


p.467_476 


p.43_61 


p.F148Y 


p.H194Q 


p.G310S 


change 


fsX459 


fsX459 


fsX613 


fsX608 


duplOA 


del19 








NLS (exon 9C) 








+' 


+' 


+ 


+ 


+ 


+ 


TAD (exon 9C) 




Involved 


Involved 


Involved 


Involved 


+ 


+ 


+ 


+ 


Function 


LOF 


N.E. 


N.E. 


LOF 


LOF 


Reduced 


GOF m 


GOF m 


GOF m 


Patient number 


3 


5 


1 


3 J 


2 


1 


1 


2 


1 558 


Occurrence 


Familial 


Familial 


Sporadic 9 


Familial 


Sporadic 


Sporadic 9 


Sporadic 


Familial 


Sporadic 


Facial features 13 


3/3 


5/5 


+ 


3/3 


0/2 




+ 


2/2 


+ 100% 


Nasal voice c 


2/3 


5/5 


N.D. 


3/3 


0/2 




+ 


0/2 


+ 32% 


Cardiovascular 


1/3 


0/5 


+ 


2/3 


2/2 


+ 


+ 


0/2 


+ 57% 


anomalies 


Hypopara- 


1/3 


3/5 


+ 


N.D. 


0/2 






0/2 


+ 60% 


thyroidism d 


Hypoplastic 


1/2 e 


0/1 ' 


+ 


N.D. 


0/2 






N.E. 


+ ? 


thymus 


Susceptible to 


N.D. 


0/5 


N.D. 


N.D. 


0/2 




N.D. 


N.D. 


N.D. ? 


infection 


Developmental 


0/3 


5/5 


N.D 


0/3 


0/2 






1/2 


38% 


retardation 


Reference 


2 


This study 3, 4 


5 


4, 6 


7 


2 


8 


2 1 


In addition to the mutations listed in this table, several missense variants and in-frame indels with unknown functions have been found in patients with isolated 



cardiovascular anomalies and in those with DGS/VCFS-like phenotype [4]. 
NLS: nuclear localization signal; TAD: transactivation domain; LOF: loss-of-function; N.D.: not described; N.E.: not examined; GOF: gain-of-function; Del: deletion; and 
Dup: duplication. 
'According to NMJD80647. 
b Suggestive of 22q11.2 deletion syndrome. 
c Velopharygeal insufficiency. 
d Hypocalcemia is included. 

e Two of the three subjects have been examined for hypoplastic thymus. 

f One of the five subjects has been examined for hypoplastic thymus. 

9 These two mutations have been inherited from apparently normal mothers. 

'The c.1 293-1 31 5del23 has been described as c.1 320-1 342del23 in the original report [5]. 

'Although the natural NLS has been disrupted, a new NLS-compatible motif (RGRRRRCR) has been created on the added amino acid sequence. 
'Another deceased individual in this family also has similar clinical features. 

k These two mutations have been identified in TBX1 analyses for patients with cardiovascular anomalies only. 
'The mutant protein is aggregated in the cytoplasm and the nucleus. 
m Gain-of-function effects have been found by in vitro studies [8]. 
doi:1 0.1 371 /journal.pone.0091 598.t002 



phenotype. Application of such powerful methods will serve to 
identify a causative gene in genetically heterogeneous disorders. 

Supporting Information 

Figure SI Comparison of amino acid sequence of 
human TBX1C and mouse Tbxl. The T-box is highlighted 
in yellow, and the nuclear localization signal in light blue. The 
region for transactivation domain is surrounded by squares. The Y 
highlighted in red denotes the amino acid residue where the 
frameshift mutation in this family has taken place. 
(TIF) 

Figure S2 Analysis of c.2968A>G SNP (p.R990G, 
rsl042636) with a gain-of-function effect in exon 7 of 



CASR. The SNP pattern is not co-segregated with the presence or 

absence of hypocalcemia. 

(TIF) 

Table SI Summary of heterozygous non-synonymous 
variants. 

(PDF) 

Table S2 A list of variants that are present in groups 
1+2 and absent from group 3. 

(PDF) 

Table S3 A list of variants that are present in group 1 
and absent from groups 2+3. 

(PDF) 
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Table S4 A list of variants that are present in group 1 
and absent from group 2. 

(PDF) 

Table S5 A list of variants that are present in group 2 
and absent from group 1. 

(PDF) 
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